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Abstract 

We study summation of sequences and integration in the quantum 
model of computation. We develop quantum algorithms for computing 
the mean of sequences which satisfy a p-summability condition and for 
integration of functions from Lebesgue spaces Lp ([0, l]'^) and analyze 
their convergence rates. We also prove lower bounds which show that 
the proposed algorithms are, in many cases, optimal within the set- 
ting of quantum computing. This extends recent results of Brassard, 
H0yer, Mosca, and Tapp (2000) on computing the mean for bounded 
sequences and complements results of Novak (2001) on integration of 
functions from Holder classes. 



1 Introduction 

Quantum algorithms and complexity are by now well studied for various dis- 
crete problems. This includes such milestones as Shor's (1994) factorization 
and Grover's (1996) search algorithm. Much less is understood about nu- 
merical problems, computational problems of analysis. These problems are 
typically defined on a continuum and/or take values in a continuum, such 
as the field of real or complex numbers, domains in finite dimensional vector 
spaces or even infinite dimensional normed spaces like function spaces. 

First results related to this direction concern the counting problem 
(Boyer, Brassard, H0yer, and Tapp, 1998) and the computation of the mean 
(Grover, 1998, Brassard, H0yer, Mosca, and Tapp, 2000) of finite sequences 
which satisfy a uniform bound (e.g. whose elements belong to the interval 



[0, 1]). Matching lower bounds were obtained by Nayak and Wu (1999) us- 
ing the polynomial method of Beals, Buhrman, Cleve, and Mosca (1998). 
Abrams and Williams (1999) proposed certain quantum algorithms for in- 
tegration. Novak (2001) was the first to provide quantum integration algo- 
rithms with matching upper and lower bounds. He studied an important 
class of integrands - functions which belong to Holder spaces. His work 
is closely related to information-based complexity - a frame in which the 
complexity of numerical problems is studied (in the classical setting). 

Mainly due to efforts within this theory, by now for many important 
problems of numerical analysis matching upper and lower complexity bounds 
(or in other words, optimal convergence rates) are known for both the clas- 
sical deterministic and randomized setting. It is a challenging task to study 
these problems in the setting of quantum computation. Once such results 
are obtained, one can compare them to the deterministic and randomized 
classical ones to understand the possible speedups by quantum algorithms. 
Novak (2001) did the first step toward this, and the present paper as well as 
related work, Heinrich and Novak (2001a,b) and Heinrich (2001), go further 
along this line. 

In the present paper we study quantum summation of sequences satisfy- 
ing p-summability conditions. These classes are larger than that of uniformly 
bounded sequences (precise definitions are given in section 3) and cannot be 
handled by the previous algorithms. But the solution of this problem is 
needed for the understanding of quantum integration in various function 
spaces (different from Holder classes) characterized by p-integrability con- 
ditions, such as the Lebesgue spaces Lp{[0, 1]*^), studied here in section 5, 
and the Sobolev spaces analyzed in Heinrich (2001). In the present paper 
we therefore develop quantum algorithms for computing the sum of such 
sequences. We also prove lower bounds which are, in many cases, matching 
with the obtained upper bounds, showing the optimality of the algorithms. 
(The picture is completed in Heinrich and Novak 2001b, where the case 
is settled which is left open here.) These results enable us to completely 
determine (in one case, up to a logarithmic factor) the optimal order of 
convergence of quantum integration in Lebesgue spaces Lp{[0, 1]'^). 

Comparing the result both for summation and integration with the ran- 
domized classical setting, we observe a considerable gain by quantum com- 
puting - the quantum speed of convergence equals the square of the ran- 
domized classical one. The gain over deterministic classical algorithms can 
even be exponential (see the details in sections 5 and 6). 

To put the problem formulations and the results on a firm mathematical 
basis it was necessary to extend the usual model of quantum computation 
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(we follow Beals, Buhrman, Cleve, and Mosca, 1998) to the setting of nu- 
merical problems, to the fields of real or complex numbers, normed spaces 
of functions etc. This extension was widely inspired by the approach of 
information-based complexity theory to numerical problems in the classical 
settings and can be viewed, in fact, as a quantum setting of this theory. 

The paper is organized as follows. The general approach is presented 
in section 2. Upper bounds for summation of p-summable sequences and 
respective algorithms are contained in section 3. General results concerning 
lower bounds as well as their application to summation are given in section 4. 
Section 5 is devoted to the application of the previous results to integration 
of functions from the Lebesgue spaces Lp{[0, l]'^). Finally, in section 6 we 
give comparisons to results in the classical deterministic and randomized 
settings and comment on some further related issues. 

For background reading in quantum computing we refer to the surveys 
Ekert, Hayden, and Inamori (2000), Shor (2000), and the monographs Pit- 
tenger (1999), Gruska (1999) and Nielsen and Chuang (2000). For notions 
and results in information-based complexity theory see the monographs 
Traub, Wasilkowski, and Wozniakowski (1988), Novak (1988), and the sur- 
vey of the randomized setting Heinrich (1993). 

2 A General Framework for Numerical Quantum 
Algorithms 

We are given nonempty sets D, K, a nonempty set F of functions on D 
with values in K and a function S from _F to a normed space G. By a 
normed space we always mean a normed linear space over K, where K is 
either R or C, the field of real or complex numbers. We seek to compute 
(approximately) S{f) for f G F, where / can only be accessed through its 
values (that is, we assume that / is given as a black box - given t E D, this 
black box returns f{t) € K). 

This general framework includes, on one hand, the binary case, where 
D = {0, . . . , N — 1}, K = {0,1}, F consists of all Boolean functions, i. e. all 
functions from D to K, and S maps F to G = R (which contains {0, 1}). 
On the other hand, in numerical problems, D is usually some subset of 
R*^, K = K, F is usually a subset of a normed linear space of functions 
(or tuples of functions) from I? to K, and 5 is a mapping (also called the 
solution operator) from i*" to G, where G is either K or a normed space of 
functions. 

We want to study algorithms and complexity of solving these problems 
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on a quantum computer. For this purpose, we adopt standard notation of 
quantum computing. Let Hi be the two-dimensional complex Hilbert space 
C^, {eo,ei} its unit vector basis, let 

be the Hilbertian tensor product of m copies of Hi. We use the standard 
identifications such as writing or \i) for e^o (8) • ■ • <8) ej^_j, where i = 
YlT=o 3><:^^~^~^ is the binary expansion of i. When identifying H^ with 
Hmi ® ■ ■ '®Hmi , where mj = m, we also identify Cj with the respective 

Cj^ (8) • • • <8) and . . . \i£), and finally also i itself with (ii, . . . , ii) in the 
respective way. For convenience we use the following notation: 

Z[0,N) := {0,...,iV-l} 

for iV G N (as usual, we let N = 1, 2, . . . , Nq = N U {0}). Let Cm = : 
i £ Z[0, 2™)} be the set of basis vectors of H„i, also called classical states, 
or basis states, and let U{Hm,) denote the set of unitary operators on Hm- 

First we introduce the notion of a quantum query (in our setting of 
D, K, F, G and S). A quantum query on F is given by a tuple 

Q = {m,m,m",Z,T,(3), (1) 

where m,m',m" G N, m' + m" < m,Z C Z[0, 2"*') is a nonempty subset, 
and 

(3 : ^ Z[0,2"^") 

are arbitrary mappings. The meaning of these components will be explained 
below. Such a tuple Q defines a query mapping (we use the same symbol 

Q) 

Q:F^ U{Hm) 

f^Qf 

as follows: Let any h £ Cm be represented as h = \i)\x)\y) with \i) G 
Cm', \x) G Cm", \y) G C-m-rn'-rn" (if m = m' + m" , we drop the last compo- 
nent). Then Qf is the unitary operator defined uniquely by its action on 
C ■ 
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where here and in the sequel always means addition modulo the respective 
power of 2, here modulo 2"^ . Let m{Q) denote the first component of Q, 
that is, the total number of qubits. If m{Q) = m, we also say that Q is an 
m-qubit quantum query. 

This notion contains the binary black box query typically used in quan- 
tum computation (see, e.g. Beals, Buhrman, Cleve, and Mosca, 1998) 
as a particular case: Such a binary query associates to a {0, l}-valued 
function k defined on Z[0, 2"* ) the operator which maps \i) \x) \y) to 
\i) \x /^(^)) \y), where i G Cm' and x E Ci. In our situation we have to deal 
with two more general domains: D and K. The mapping r : z — ^ T{i) G D 
describes the (chosen by the algorithm designer) correspondence of binary 
strings with certain elements of the domain of definition of functions from 
F. Since at request r{i) the black box returns f{T{i)), which is an element 
of K, we need a second mapping P, which maps ("codes") elements from 
K into binary strings. (This is also chosen by the algorithm designer.) As 
usual, the untouched part \y) stands for "working bits". 

Note that, by the definition, a quantum query on F is also a quantum 
query on any other nonempty subset Fi C T[D,K), and in particular, on 
J-{D, K) itself. Here J-{D, K) denotes the set of all functions from D to K. 
Indeed, the mapping Q/ is defined for each / G J^{D, K). 

Next we define quantum algorithms in the general framework of D, K, F, 
G and S. It will be convenient for us to introduce algorithms with multiple 
measurements. We show later in this section how they can be simulated 
by algorithms with one measurement. Let us first describe informally what 
we mean by a quantum algorithm with k measurements: Such an algorithm 
starts with a fixed basis state 6o and applies in an alternating way unitary 
transformations (not depending on /) and a certain query, associated to the 
algorithm. After a fixed number of steps the resulting state is measured, 
which gives a (random) basis state ^o- This state is memorized and then 
transformed (e.g. by a classical computer) into a new basis state 6i. This 
is the starting state to which the next sequence of quantum operations is 
applied (with possibly another query and number of qubits). The resulting 
state is again measured, which gives the (random) basis state ^i. This state 
is memorized, and 62 is computed from ^0 and ^1, and so on. After k such 
cycles, we obtained • • • ) ^fc-i- Then finally an element of G - the output 
of the algorithm - is computed (e.g. again on a classical computer) from the 
results of all measurements: ip{^o, ■ ■ ■ , Cfc-i)- 

Now we formalize this: A quantum algorithm on F with no measurement 
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is a tuple 

A = {Q,{U,)^^,), 

where Q is a quantum query on F, n € Nq and Uj € lA{H„i) (j = 0, . . . , n), 
with m = m{Q) (in the case n = 0, no query Q is needed). Given such an 
A and / G F, we let ^dj G U{Hm) be defined as 

Af = UnQfUn-l...UiQfUo. (3) 

We denote by ng(^) := n the number of queries and by m{A) = m = m{Q) 
the number of qubits used by A. We also introduce the following notation. 
Let Af{x,y) for x,y G Z[0,2'") be given by 

A}\y)= M^^y)\^)- (4) 

xez[o,2'") 

Hence {Af{x,y))x,y is the matrix of the transformation ^/ in the canonical 
basis Cm- 

A quantum algorithm on F with output in G (or shortly, from F to G) 
with A; measurements is a tuple 

A={{A,)tl{bi)tl^)^ 

where /c G N, and (^ = 0, . . . , A; — 1) are quantum algorithms on F with 
no measurements. To explain the other components, set mi = m{A(). Then 

6oGZ[0,2"*o), 

for 1 < ^ < A; — 1, 6^ is a function 

i=0 

and is a function with values in G 

k-l 

]JZ[0,2™*)^G. 

e=o 

We also say that A is a quantum algorithm with measurement (s), or just a 
quantum algorithm. 

Let Vo{G) denote the set of all probability measures on G whose support 
is a finite set. The output of A at input f & F will be an element A{f) G 
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Vo{G) (we use the same symbol A for the mapping A : F ^ Vq{G)). We 
define A{f) via a sequence of random variables i£,ij)£Zo assume that 
all random variables are defined over a fixed - suitably large - probability 
space (il, S, P)). So let / S F be fixed. Now let ^ij be such that 

P{Co,f = x} = \Aoj{x,bo)\^ (5) 

and, for 1 < ^ < A: - 1, 

P{C£,/ = x\^oj = xo,.. . ,6-1,/ = Xi-i} = \Aij{x,bi{xo,.. . (6) 

Clearly, this defines the distribution of (6j)^~q uniquely. Let us define for 
xoGZ[0,2™«),...,Xfc-ieZ[0,2-'=-i) 

PAj{xo,. . . ,Xk-i) = \Aoj{xo,bo)\'^\Aij{xi,bi{xo))\'^ . . . 

. . . |Afc_ij(xfc_i,6fc_i(xo, . . . ,Xfc_2))P- (7) 

It follows from (|5|) and (|6|) that 

P{6,/ = xo, . . . , 6-1,/ = Xk-i} = PAj{xo, Xk-i). (8) 

Finally we define the output of A at input / as 

^(/) = dist((^(eo,/,...,6-i,/)), 

the distribution of (/?(,^o, /)•••) 6-i,/)- This random variable takes only 
finitely many values in G, hence the support of A{f) is finite (and no mea- 
surability problems related to the target space G will arise). It follows from 
(^) that for any subset C C G 

Mf){C}= PAj{xo,--.,Xk-i). (9) 

We note that, analogously to quantum queries, a quantum algorithm on F 
is automatically also a quantum algorithm on any nonempty Fi C J-'{D, K). 

The number nq{A) := J2e=o ^gi^e) is called the number of queries used 
by A. This is the crucial quantity for the purposes of our query complexity 
analysis. (In section 6 we give some comments on the cost in the bit-model.) 

Let < 9 < 1. For an algorithm A as above we define the (probabilistic) 
error at / G -F as follows. Let ^ be a random variable with distribution 
A{f). Then 

e{S,A,f,e) = inf {e I P{\\S{f) - C\\ > e} < 6} 
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(note that this infimum is always attained). Hence e{S,A,f,9) < e iff the 
algorithm A computes S{f) with error at most e and probability at least 
1-61. We put 

e{S,A,F,9) = sup e{S, A, f, 9) 

(we allow the value +00 for this quantity). Furthermore, we set 

e{S,A,f)=e{S,A,f,l/4) 

and similarly, 

e{S,A,F)=e{S,A,F,l/A). 

The central quantity of our study is the ra-th minimal (query) error, defined 
for n G No by 

e^{S,F) = mf{e{S,A,F) \ A is any quantum algorithm with nq{A) < n}, 

that is, the smallest error which can be reached using at most n queries. 
The query complexity is defined for e > by 

compf(5,F) = 

mm{nq{A) \ A is any quantum algorithm with e{S,A,F) < e} 

(we put comp|(5, F) = +00 if there is no such algorithm). It is easily 
checked that these functions are inverse to each other in the following sense: 
For all ra G No and £ > 0, en{S,F) < e if and only if comply (5*,^) < n 
for all ei > s. Hence it suffices to determine one of them. We shall usually 
choose the first one. 

Our first general result shows the tight relation between algorithms with 
several measurements and (the conceptually simpler) algorithms with one 
measurement. It states that an algorithm with several measurements can 
always be represented equivalently by an algorithm with one measurement 
and twice the number of queries (at the expense of an increased number of 
qubits) . 

Lemma 1. For each quantum algorithm A from F to G with k measure- 
ments there is a quantum algorithm A from F to G with one measurement 
such that nq{A) = 2nq{A) and 

m = A{f) 

for all f £ F. 
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Proof. By '£-th quantum cycle' we mean the quantum operations in the 
original algorithm before the first measurement if £ = 0, and between the 
^-th and the i + 1-st measurement if 1 < i < k — 1. The idea of the proof is 
easy: We simulate the k queries by one query and instead of intermediate 
measurements we 'store' the results of the cycles in different components 
until the final measurement (a pseudo-code is given below). Let us now 
formalize this and check that the corresponding probabilities coincide. Let 
the original algorithm be given by 

^=((Mto>(Mto.^)' 

where 
and 

Qi = (m£, m'^, m'l, Zi, r^, (3^). 

By adding, if necessary, qubits, which are set to zero and remain so during 
the whole ^-th cycle we may assume without loss of generality that m'^ = m! . 
LetQ ko = [log/c] , define m' = m' + /cq, and Z C Z[0, 2™') by 

Z = |0 < £ < A; - 1, i G ZJ. 



Now we define 



Moreover, we set 



t: Z ^ D 

for {£,i)eZ. 



k-l 

rri' 

1=0 

P : K ^ Z[0,2^") 
P{s) = {Po{s),...,Pk-i{s)) for seK. 

k-l 

in = ko + in" + 

e=Q 

Q = {in, in' , in" , Z, r, /?). 
Let us fix the following notation: Consider the splitting 

^Throughout this paper log stands for logj. 
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The representation of a basis state 

\i) \u) \xo) ■ ■ ■ \xk-i) 
refers to this spUtting. We also need refined sphttings. We represent 

Hfh" = H^ii (g) • • • (g) Hjnn_^ , 

and 

\u) = \uo) . . . \uk-l) 
corresponds to that spHtting. Similarly, 

Hrrii = H^i (g) H^u (g) 

with the respective 

= \h) \yt) W) . 

Next we define the following unitary operators on iJ^ by their action on the 
basis states: 

J |i) |n) |xo) . . . \xk-i) = \i) \Qu) \xo) . . . \xk-i) , 

where means subtraction modulo 2"^ and Qu stands for OQ u, 

C \i) \u) \xq) . . . \xk-i) = K e 1) \u) \xq) . . . \xk-i) , 

(OT £ = 0,...,k-l, j = 0,...,ne, 

Ti \i) \uo) . ■ ■ \u£) . . . \uk-i) \xo) . . . \yi) \zi) . . . \xk-\) 
= \i) \uo) . . . \ue) . . . \uk-i) \xo) . . . \ie) \ye ® ue) \ze) . . . \xk-i) 

Uej \i) \u) \xo) . . . \xe) . . . \xk-i) = \i) \u) \xo) . . . {Uej \xe)) . . . \xk-i) , 

Pg \i) \u) |xo) . . . \ye) \zi) . . . \xk-i) 
= \i) \ie) \u) \xo) . . . \ye) \ze) . . . \xk-i) , 

and finally, for £ = 1, . . . , A; — 1, 

Be \i) \u) \xo) . . ■ \xi_i) \xi) . . . \xk-i) 

= \i) \uo) \xo) ■ ■ ■ \xe-i) \x(, e bi{xQ, xe-i)) . . . \xk-i) . 
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Now we present the simulation of the queries Qtj hy Qf. Let < i <k — \. 
It is readily checked, that if we apply the operator P^^QfP^ to the state 

\i) |0) \xo) ...\xk-i) , 

we get 

\£)\p{f{T{e,ie))))\xo) \xk-i) 

= \e) mfMieW . . . \Pk-iifire{iem ko) I^it-i) , 

provided i£ & Z^. Applying then Ti to this state gives 



1^) i^)))^ \xo) . . . \ie) \ye © PiifiTiiie)))) \zi) . . . \xk-i) 

= \t) |^(/(r(£, h)))) \xo) . . . {Qej \xe)) . . . \xk-i) . 
Next J is applied which yields 

1^) \eP{fir{£,ie)))) ko) • • • (Qej \xe)) ■ ■ ■ \xk-i) , 

and finally the appHcation of Pf^QfPe produces 

1^) |0) |xo) . . . (Q^j Ix^)) . . . |a;fe_i) . 

If i£ ^ Z£, this also holds, which is checked in the same way. Hence we 
showed that 

Pf^QfPiJTePf^QfPe \£) |0) \xo) . . . \xi) . . . \xk-i) 

= \£)\0)\xo)...{Qi,f\xi))...\xk-i). (10) 

The new algorithm can now be described as follows: 

initialize |0) |0) |5o) |0) . . . |0) 
for^ = 0,...^A;-l do 

apply U^^Q (beginning of £-th. cycle of original algorithm) 

for j = 1, . . . ^ 

apply P^^QfPeJTePf^QfPe 

apply U£j (end of £-th cycle of original algorithm) 
if ^ 7^ A; - 1 

apply -B^+i (computing as initial state of next 
cycle) 
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apply C (increasing the counter by one) 
measure all qubits corresponding to the components 
Hmo, ■ ■ ■ ,Hmk-i (let |xo) . . . \xk-i) be the result) 
compute (p{xo, . . . , x^-i)- 

The starting passage through the outer loop {£ = 0) acts as follows: 
|0) |0) |6o) |0)('-^) ^ |1) |0) (^Y^AojixoM \xo) Mxo))^ \of-'^ • 

The passage with index i, 1 < i < k — 2, maps each basis state of the form 
|^)|0)|xo)...|x,_i)|y)|0)('=-^-^) 

into 



\e + 1) |0) |xo) . . . \xe-i) I Aejixe, y) \xi) \b. 



,+i(xo,...,x,)) |0)('=-^-2) 



Finally, the last passage {(. = k — 1) acts as follows: 
|A;-l)|0)|xo)...|xfc_2)|y) ^ 

\k - 1) |0) |xo) . . . \xk-2) ^ Ak-ij{xk-i,y) \xk-i) 

From this it follows that the overall result of the algorithm before measure- 
ment is the state 

YZ ^o,/(a;o,&o)-4ij(xi,6i(a;o)) . . . 

. . . j(xfc_i, 6fc_i(xo, . . . , Xk-2)) |0) |xo) . . . \xk-i) . 

The probability of measuring |xo) • • • is thus 

l^o,/(a:^o, bo)\^\Aij{xi,bi{xo))\^ . . . \Ak-ij{xk-i,bk-i{xo, . . . , Xk-2))\^, 
which equals 

P{?0,/ = Xq,..., ik-lj = Xk-l], 

by (§) and (|). This proves the lemma. □ 
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We will sometimes write that we repeat a quantum algorithm a number 
of times, or, more generally, that we apply to / € _F a finite sequence of 
algorithms Ai from F to Gi (i = 0, . . . , M — 1) and combine the results by 
the help of a classical computation. Let 

ip : Go X ■ ■ ■ X Gm-1 G 

be any mapping. Using our notion of a quantum algorithm with measure- 
ments, a formal representation of the composed algorithm A, which we write 
symbolically as 

A = V(^o,---,^M-i), (11) 
can easily be given as follows: Let 

put k = X^^o ^ ki, let the set 

T = {{i,t) I i = 0,...,M-l, ^ = 0,...,A;i-l} 
be equipped with the lexicographical order, and let 

ip = Ipiifo, . . .,(pM-l)- 

Then we define 

il^{Ao, . . . , Am-i) = {{Ay)ver, {by)y^r,^)- 

The next lemma gives some further description of the composition and 

is readily checked using the definition of a quantum algorithm. We need 
the following notation: For probability measures //q, • • • , A*M-i £ 'PoiG) let 
V'(y^O) • • • ) A*M-i) £ 'PoiG) be the measure induced hy /iqX ■ ■ ■ x /im-i via 
on G, that is, for C C G, 

i){no,...,liM-i){C) = {noX-- - X HM-i){'4^'^{C)). 

Lemma 2. For each f & F, 

i;{Ao, . . .,AM-i){f) = i^iAoif), Am^iU)), 

or stated equivalently, if iCi)iLo^ independent random variables with dis- 
tribution Ai{f) respectively, then 

V'(Ao, • • • , ^M-i)(/) = dist(V'(Co, ■ • • , Cm-i)). 
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Moreover, 

M-l 

nq{ip{Ao, Am-i)) = ^ nq{Ai). 

1=0 

The next lemma concerns the special case of repeating an algorithm. 
It describes a standard technique of boosting the success probability. For 
completeness, we include the short proof. Let G = R, M G N and denote by 
'00 : —>■ R the mapping given by the median, that is, ipo{C'0^ ■ ■ ■ ■> c^M-i) 
is the value of the of the \{M + l)/2]-th element of the non-decreasing 
rearrangement of (oj). For any algorithm A from F to R denote tpQ{A^) := 

Lemma 3. Let A be any quantum algorithm and S be any mapping from F 
to R. Then for each f E F, 

eiS,MA^)J,e-''/^)<<S,A,f). 

Proof. Fix f E F. Let Co) • • • ; Cm-i be independent random variables with 
distribution A{f). Let Xi be the indicator function of the set {\S{f) — d] > 
e{S,A,f)}. Then P{xi = 1} < 1/4. Hoeffding's inequality, see e.g. Pollard 
(1984), p. 191, yields 

{M-l ^ fM-l -\ 

> M/2\ < P I E - EXi) > M/a\ < e-^l\ 

Hence, with probability at least 1 — e"-^/®, 

|{i||-S(/)-a|<e(5,A/)}|>M/2, 

which implies 

l'5(/)-Vo(Co,...,CM-i)|<e(5,A/)- 

□ 

Another way of building new algorithms from previous ones will also be 
important Jor us^To explain it, let / F C T(p, /C) and / F C Tip, K), 
where D, D, K, K are nonempty sets. In the construction of a new algorithm 
A on F we sometimes construct from / a function / = T(J) £ F to which 
we want to apply an already developed algorithm A on F. By definition, 
the algorithm A on. F can only use queries Q on F itself, while we need to 
use Qr(/)) where Q is a query on F. Nevertheless often a solution can be 
found as follows: We simulate Qr{f) either as Qf with a suitable query Q 
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on F or as Bf, where B is an algorithm without measurement on F. The 
details are given below. 

The first result covers the simple situation where one query is just re- 
placed by another. Let t] : D ^ D and g : K ^ K he arbitrary mappings 
and define T : F ^ F hy 

^(/) = ^^o/or?. (12) 

Lemma 4. Let F be a mapping of the form (fJ^j. Then for each query Q 
on F there is a query Q on F such that m{Q) = m{Q) and for all f ^ F 

Qf = Qr{f)- 

Proof. Let 

Q = {fh, m', m", r, (3). 

Then we define 

Q = (m, m', m", Z, r, 

where t = r] or and (3 = (3 o g. Now the lemma follows directly from the 
query definition. □ 

The second result in this direction is slightly more technical. We assume 
that we are given a mapping F : — > i*" of the following type: There are an 
m* G N and mappings 

ri : D 

(3 : ir^Z[0,2™*) 

g : 5xZ[0,2™*)^^ 

such that for f G F and s £ D 

nf)is) = g{s,/3ofor^{s)). (13) 

Lemma 5. Let Q be a quantum query on F and let T be a mapping of the 
above form (jl^. Then there is a quantum algorithm without measurement 
B on F such that nq{B) = 2, m{B) = m(Q) + m* and for all f £ F, 
X G Z[0,2"*(<3)), 

i?/|x)|0)„, = (Qr(/)k))|OL*, 
where |0)^, stands for the zero state in Z[0, 2™ ). 
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Proof. Let _ 

Q = {m,m' ,m'' , Z,t,P), 



and put 



m = m + m*, m! = m' , m" = m* , 



Z = Z, T = rj OT, 
let P he as above and define 

Q = (m, m', m" , Z, r, /?). 

We represent 

Hm = Hfj^, (g) Hffi,/ (g) Hffi_fj^,_ffi,/ (g) Hm* , 

a basis state of which wih be written as 

\i) \x) \y) \z) . 
Define the permutation operator P by 

P\i) \x) \y) \z) = \i) \z) \x) \y) , 
the operator of sign inversion 

J\i)\z)\x)\y) = \i)\Qz) \x)\y), 

and finally 

T \i) \z) \x) \y) = \i) \z) x®j5o Q{T{i), z)^ \y) 

a i & Z, and 

T\i)\z)\x)\y) = \i)\z)\x) \y) 
ii i ^ Z. We define B by setting for f £ F, 

Bf = p-^QfJTQfP. 

Let us trace the action of Bf on 

\i) \x) \y) |0) . 

First we assume i e Z. The transformation QfP leads to 

\i) Wifirim k) \y) = l'^) 1/3 o / o r? o |x) \y) 
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Then the above is mapped by T to 

\i) \I3 o f o r] o T{i)) X ® (5 o giji^i), f3 o f o rj o t(«)) 

and P^^QfJ gives 

\i) \x e mfWm) \y) |0) = (Qr(/) N> \x) \y) ) |0) . 
The case i ^ Z is checked analogously. □ 



Corollary 1. Given a mapping T : F ^ F as in ( |7^ j or JT^j, a normed 
space G and a quantum algorithm A from F to G, there is a quantum algo- 
rithm A from F to G with 



nq{A) 

and for all f ^ F 



I nq{A) in case of ( f7^ 
I 2nq{A) in case of (jil 

A{f) = A{T{f)). 



Consequently, if S : F ^ G is any mapping and S = S o T, then for each 
n e No 



el{S,F) < el{S,F) in case of g^, and 
eUS,F) < el{S,F) in case of 1^ 



Proof. Let 

and fhe = m{Ai). Then fov f e F, < £ < k, 



In case of (|12| ) we obtain A by just replacing Qg by Qg from Lemma It 
follows from (0) and @ that 

A{f) = I(r(/)). 
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In case of ( p^ ) we replace Qn by from Lemma |5|, C/^j- by f/^j = C/^j (81 
Mh^,, where Idn^, is the identity on Hm*, the state 60^ by |6o) = 

60^ |0)m* ^^'^j for 1 < ^ < — 1, the mappings 

be : J|Z[0,2'^0 ^ Z[0,2^0 

i=0 

by 

e-i 

be : Yl (Z[0,2™0 X Z[0,2™*)) ^ Z[0,2'^0 x Z[0,2™'), 

defined by 

be{{xo,yo), . . . , (xf_i,y£_i)) = (^^(xo, . . . ,x<?_i),0). 
Finally, we replace 

k-l 

if: J|Z[0,2"^) 

by 

k-l 

(Z[0,2"^) X Z[0,2"*)) -^G, 

defined as 

^iixo,yo), {xk-i,yk-i)) = 'fixQ, ■ ■ ■,xk-i). 

It follows that 

A f((x v) (z 0)) = ( if y = 

Aej[(x,y),[z,U)) 1^ otherwise, 

and therefore, by (|^, 

PAji{xo,yo) {xk-i,yk-i)) 

^'A,r(/)(^o, • • • if yo = • • • = yfc-i = 



otherwise, 
which together with (^) yields 

A{f) = A{T{f)). 

This proves the first part of the statement. The second part is an obvious 
consequence. □ 
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Finally we state some elementary but useful properties of e^- For A G K 
define XS : F ^ G hy (A5)(/) = XS{f) {f £ F). Furthermore, in the case 
K = Kwe denote XF = {Xf\f e F}. 

Lemma 6. Let S,T : F ^ G be any mappings, n G Nq and assume that 

en{S,F) is finite. Then the following hold: 

(^) 

el{T,F) < eUS,F) + snpWnf) - Sif)\\. 

f<^F 

(a) For each A € K 

el{XS,F) = \X\el{S,F). 

(Hi) If K = and S is a linear operator from J-{D, K) to G, then for all 
A G K 

el{S,XF) = \X\el{S,F). 

Proof. The first two statements are simple consequences of the definitions. 
Let us verify the third one. Let F = XF, T : F ^ F he defined as r(/) = A/, 
which is of the form (12). We assume A 7^ 0, the case A = follows trivially 
from (ii). Since S is linear, we have 

X-^SoT = S, 

and hence, by Corollary |^ and statement (ii) above, 

el{S,F) < \Xr'eUS,F) = \X\-'eUS, XF). 

Replacing F by XF and A by A~^, we get 

el{S,XF)<\X\el{S,F), 

which completes the proof. □ 



3 Quantum Summation 

In this section we study summation of sequences or, what is essentially the 
same, the computation of the mean, on a quantum computer. For a fixed 
iV G N we set L> = Z[0, iV), K = R, G = R, and for 1 < p < 00 let 
denote the space of all functions / : D — > R, equipped with the norm 

/ N-l 
\ i=0 



19 



a p < oo and 

(Note that Lp is just the space Lp{D, ij,), where is the equidistribution on 



TN = max \f(i) 
^'^ 0<i<N-l ' ^ ' 



D.) Define 5jv : ^ R by 



N-l 



i=0 



We let 



F = := B{L^) = {feLi 



<1} 



be the unit bah of Lp . We also define 



B. 



N 



{f : D^K\0< f{i) < 1 for ah i} 



and 



^^,0 = {/ : - {0, 1}} . 
When we consider B^ we put K = {0, 1}. Clearly, 



B, 



N 



00,0 ^oo,+ ^1} 



2N ^ t2N 



whenever 1 < q < p < oo. Therefore, we will also consider Sn as acting on 
"^00,+ • We use the following standard representations depending 
on the range of / G -P: Given a, 6 G R, a < 6, and k G N, define (3K,a,b '■ 
R^Z[0,2'^)by 



l3K,a,b{x) = < 
So fov a < X < b 



( 2'^ -I if X > 6 

if a; < a 

1 if Iff G[^,^),zGZ[0,2«). 



(14) 



PK,a,b{x) 



2K 



X — a 
b — a 



and hence, ioi a < x < b, 

a+{b-a) 2-^(3^,a,b{x) <x<a+{b-a) 2-'^(/3«,„,6(x) + 1). (15) 

First we state the basic result on quantum counting due to Brassard, H0yer, 
Mosca, and Tapp (2000). 
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Lemma 7. There is a constant c > such that for all n, iV G N there is 
a quantum algorithm A from B^q to R such that nq{A) < n and for each 



oo.O 



Remark. Throughout this paper we often use the same symbol for possibly 

different constants. These constants are either absolute or may depend 
only on p - the summability parameter of the Lp-spaces considered (in all 
lemmas and theorems this is precisely described anyway by the order of the 
quantifiers) . 

Proof. We refer to Brassard, H0yer, Mosca, and Tapp (2000) for details 
of the algorithm, its analysis and the resulting estimates. For us, there 
remains one detail to be verified. Their algorithm makes use of the controlled 
application of the Grover iterate and assumes that an implementation of 
this procedure is available. This means, roughly, if Y stands for the Grover 
iterate, we must be able to implement an operation which maps an element 
\i) \ k) to (Y^ \ \ k) (that is, different basis elements may be subject to 
different powers of Y). Since Y involves a query call, it is not immediately 
clear, how this could be achieved within the rules delveloped in section 2, 
that is, in our model of computation and its way to use queries. So we 
supply the needed argument here. It is a simulation procedure, similar to 
the ones above. 

The parameters of the algorithm will be the following. It has one mea- 
surement, and the query Q is determined by 

m' = [log N~\ , m" = 1, m* = [log n] , 
m = m' + 2m* + 2, Z = Z[0,N), 
t:Z^ Z[0, T^') and /3 : {0, 1} ^ {0, 1} the identities 

(recall that K = {0, 1}). Let 
and let the basis state 

\i) \x) \j) \k) 

correspond to this splitting. Let ^n,m* be the n-term quantum Fourier 
transform on m* qubits, 

\k) otherwise. 
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Define $ G U{Hm) by 

<P\i) \x) \j)\k) = \i)\x) \j) {<Pn,m*\k)). 

Furthermore, let Vq G U{Hm') be the Walsh-Hadamard transform VFjv, if N 
is a power of 2, and let Vq = ^N,m': if not. Define Xq € U{H^i) by 



and unitary transforms on iJ^ by 



|i) if i = 
|i) otherwise, 



TK)|x)|j)|fe) 



V\i)\x)\j)\k) = {Vo\i))\x)\j)\k), 

(Xo |a;) |j) if j<k 
\i) \x) \j) \k) otherwise, 

i-ir+^\i)\x)\j)\k) if j<k 

\i) \x) \j) \k) otherwise, 
C\i)\x)\j)\k) = \i)\x)\j®l)\k). 



Now we define the algorithm as follows. For / G q set 

Yf = CVXV-^QfTQf. 
The unitary transform of the algorithm is given by 

The initial state is 

6=|0)|0)|0)|0). 

Let us now follow the action of the algorithm. The element 5 is transformed 
by into 

(yo|0))|0) |0) i^n,m* |0)). 

Note that this vector is a linear combination of basis states of the form 

|i) |0) |0) \k) 

with i < N and k < n. Next consider the application of Yf to a basis state 
of the form 

\i) |0) \j) \k) (16) 
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with i < N and k < n. First we assume j < k. Then QfTQf produces 

\i) |0) \j) \k) . 
After the appUcation of CVXV^^ we get 

{-i)m+i^v,XoV^-^\i)m Uei)|fc), 

which is a hnear combination of vectors of the form 

|i'>|0) \3 + l)\k) 
with i' < N . If j > k, the apphcation of 1/ to (jl^) gives 

|i)|o) |jei)|/c). 

It is now clear that Yf = CVXV^^Q jTQ j reahzes the Grover iterate on the 
first component if j < k, that Yp^^ is the controhed (by k) apphcation of 
it and the whole algorithm, considered just on the first and last component 
\i) \k), is the algorithm "Est_Amp" of Brassard, Il0yer, Mosca, and Tapp 
(2000), if we define f on the measured state 

\y) = \i) k) \j) \k) 

as 

(p{y) = sin^ I TT- 
\ n 

The required estimate (with a concrete value of the constant) is contained 
in Theorem 12 of that paper. Since our implementation requires 2n queries, 
we rescale n and modify the constant appropriately. □ 

The next result is essentially a translation of Lemma |^ into the setting of 
The idea of using comparison queries is due to Abrams and Williams 

(1999). 

Lemma 8. There is a constant c > such that for all VjU, N £ IS! there is 
a quantum algorithm A from to H such that nq{A) < un and for each 

f e C+ 

e{SN, A, /, 2-") < c i^/S^fn-^ + n'^^ 
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Proof. Let k G N be such that 2" > and put A''o_= N2'^. We shall apply 
Corollary || with F = _^_ and F = B^q. Let A be any algorithm from 

B^\ to R with one measurement, which satisfies the conclusion of Lemma 

^ with nq{A) ■.= n<n. Let A be given by 

I=(Io,6,^), Io = (Q,(f/,)?=o), 

with 

Q = {m,m ,m' , Z,T,i3), 

where Z C Z[0,2'^'), t : Z ^ Z[0,iVo), and (3 : {0,1} ^ Z[0,2""). We 
identify 

Z[0,iVo) = Z[0,iV) X Z[0,2'") 
and write correspondingly for z £ Z, 

T{z) = ii{z),y{z)). (17) 

Now let P = /3«,o,i as defined in (|l|). For each / E B^^+ define r(/) G 5,^" 
by setting for (i, y) G Z[0, N) x Z[0, 2"^) = Z[0, iVo 

r(/)(i,y) = 



CXD,0 



Note that 



and consequently 



1 if 2/</3(/W) 
otherwise. 



|{y:r(/)(.,y) = l}|=/3(/(i)). 



i=0 

By (H), 

5jVor(/) < 5jv/ < SNoTif) + 2-^ < SN.nf) + n-\ 

The mapping L : / ^ r(/) is easily seen to be of the form (|l3|) (with 
r]{i,y) = i and (3 as defined above). By Corollary |l| there is an algorithm A 
on + such that nq{A) = 2nq{A) and A{f) = A{T{f)). To estimate the 
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error of A, fix any / G _^_ and let C be a random variable with distribution 
^(r(/)). Then, with probability at least 3/4, 

\SNf-C\ < \SNf-SNor{f)\ + \SNor{f)-C\ 



< c I yJSNfn ^ + n 



-2 



Now we use Lemma ^ to boost the success probability by repeating A ciu 
times, where ci = [8/ log e] , and computing the median, which gives the 
desired error estimate 



e(S,A*,/,2--) < c' (yS^fn-^ + n 



for the algorithm A* = 'iI)q{A^'^^), whose number of queries is bounded by 
2civn. A scaling of n at the expense of enlarging the constant gives the 
result as required. □ 

Now we are ready to estimate the numbers et{SN,Bp). Note that this is 
nontrivial only when n < N. For n > N a classical computer suffices, or, 
to put it more formally into our framework, we have en{S]\[,Bp) = 0, since 
with N queries (and a suitable number of qubits) the sum can be determined 
up to each degree of precision by e.g. simulating a classical computation. 

The following is the main result of this section. For the sake of later 
reference we also include the already known case p = oo due to Brassard, 
H0yer, Mosca, and Tapp (2000), which we deduce formally from the case 
2 < p < oo, but which is, in fact, an immediate consequence of the previous 
two lemmas. 

Theorem 1. Let 1 < p < oo. Then there is a constant c > such that for 
aUn,N G'N, n> 2 

{n^^ for p > 2 

n~ ^ log^/^ n log log n for p=2 
n-2(i-i/p) for p < 2. 

Proof. Let 1 < p < oo. Fix k G Nq (to be specified later) and define for 
Tj = {iGZ[0,iV)||/(i)|>2'=}, 
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for cj = 0, 1, 



j;''^ = {iGZ[o,iv)|o<(-ir/w<i}, 

and for i = 1, . . . , k, 

J^^'- = |, G Z[0,N) I 2^-1 < i-irfd) < 2^} 



Note that 



hence 



Holder's inequahty together with (15) gives 



i/p' / 1 



i/p 



< 2-^"=/^''||/||^//||/||^. = 2-(^'-i)'=||/|| 



where l/p+ = 1. Furthermore, 

1 5] 2^(^-1)1^;'^ < ii/ir^ 



i<^<fe 

tT = 0,l 



which gives, in particular, 

\jjn<N2-p('~'^\\f\\l^ ii>i). 

Now define gj" G + for < ^ < A:, a = 0, 1, 



9f (0 



otherwise. 



Consequently, < g^j^ < 1, so g^j^ G Clearly, 



SNg'f < iV- V/'l 
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and 

= ^( E E /w + E/(^ 

= ^ + (23) 



N 

0<e<k ic-rk 
<T = 0,1 *^^/ 



Now the idea is to compute Sj^ig^j" by the algorithm from Lemma ^ for all 
i and 0", and from the results (in a classical way) the first sum of equation 
(p3|) . Fix f£,n£ G N (to be specified later) and let, according to Lemma 



Ai be an algorithm on -S^ _,_ such that nq{A() < uirii and for all g G ^S^^^, 
e(5iv,I^,9,2-"^) < c (x/S^n^i +n72) . (24) 
We define for x G R, o" = 0, 1, 

eoA^) - I otherwise, 
and for £ = 1, . . . ,k — 1, 

Furthermore, we let rj be the identity on Z[0, A^). Then for each f € L, 

9f = Ql,a o f or]. 



{-iy2-^x if 2^-^ < (-l)'^x < 2^ 
otherwise. 



AT 



By Corollary |l] there is an algorithm A^^a- on with 

nq{AtA = nq{Ai) 

and 

A,,„{f) = Mg'p (25) 



for all / G Lp . We define A as being composed of A^^^j (in the sense of dTT] 
as follows: 

^= E (-i)''2%,.. 



0<l<k 
o-=0,l 
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To estimate the error of A, fix any / G and let {Ce,a \ < £ < k, a = 
0, 1} be independent random variables with distribution respectively. 
Define 

C = Yl (-1)^2^0,- (26) 

o<e<k 

CT = 0,1 

It follows from Lemma ^ that 

Aif) = dist(C). (27) 
By (H) and (|2|), we have, with probability at least 1 - 2"^^ 



and therefore, with probability at least 1 — 2 J2e=o 2 



o<e<k o<e<k 

o-=0,l o-=0,l 



{-ir2\SN/f - < c J] 2^ ( ^JSN9'fnJ' + 

0<i<k 
o-=0,l 

hence, by p^) and (p6| 



i5^/-ci<c 5] 2^ +^7') + 1^ E /(^) 

a=d;i 

which gives together with (|l9|), (H) and (|2l|) 
e(5^,A,/,2j]2--0 

o<e<k ^ ' 
< c 2^(0^^^n7i + n72)+2-(^"i)^||/r,^, 



0<<<fc 



< 2c^(2(i-^/2)vi/ir^V2s-2 



(28) 



+ 2moi + 2-(^'-i)^||/||^, (29) 
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(recall the remark about constants after Lemma Moreover, we have 



(30) 



Now we choose the parameters k, vi and ni in a suitable way and prove the 
error estimates. First we consider the case 2 < p < cxd. Here we put 



k 



1 



p — 1 



log n 



Define, furthermore, = [21og(^ + 1)] + 4, hence 



i=0 



Finally, let 



2(l/2-p/4)^^ 



This together with ( |30[) implies 



(31) 



(32) 



(33) 



n, 



,(^)<2^(r21og(£+l)l+4) 



,(l/2-p/4)<! 



n 



< cm (34) 



i=0 



for some constant ci > 0. It follows from (p2|), (po]), ( pSf ) and (pi]) that 
eiSN,AJ) 

k 

< e(5iv,A/,2^2-'^0 

k 



< (2(i/2~p/4Kn-i||/||^/v' + 2P^/2n-2) + en"! + 



V 



\ V V 

< cn-imax(||/||P^,l). 



Consequently, 



6(5^,^-6^) < cn-\ 
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which together with ( |34D imphes the desired result in the case 2 < p < oo. 
Note that the case p = oo also follows since C for any p < oo. 
Now we suppose 1 < p < 2. Here we choose 



k 



p 



■ log n 



(35) 



2-(l/2-p/4)(fc-^) 



n 



(36) 



and I'i = [21og(A; — £ + !)] +4, which implies that (32) holds again. Fur- 
thermore, by (^0|), 



n, 



,(^) < 2^([21og(A;-^ + l)] +4) r2-(i/2-p/4)(fc-£) J < ^^^^ (37^ 



We get from (H), (|9|), and 



e(5;v,A/) 



(l-p/2K+(l/2-p/4)(fe-£)^-l|| ^||P/2 _^ ^e+{l-p/2)(k 



+ 2' 



"n 



+ cn^^ + 2-^f-^^'^\\f\\l^ 
p 



(l/2-p/4)(fc+£)„-l|| f ||P/2 _^ 9fc-p(fc--^)/2^-2 



^^2(i/2-p/4)fc^-i + 2-(f-i)^||/||^^ 

p 

< c (2(i-P/2)'=n-i||/||^/v' + 2'=n-2 + 2(i/2-p/4)/c^-i + 2-(p-i)'= 



p 



< cn 



'2(l-l/p) 



max 



r JV 



(38) 



Now ( ^7[ ) and (38) yield the needed result. 

Finally, we consider the case p = 2. Here we define 



(39) 



n£ = uq = |"n(logn) ^ (log log n) 
(recall that we assumed n > 2, so no is well-defined and no > 1), furthermore 

A;=riognol (40) 
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and 



It follows that 



and, by (|3g) 



vi = UQ= \\og{k + 1)1 + 3. 



^=0 



1=0 



'^g(^) < 2(A; + l)fono < cin. 



By ([4^ ) and (|2^), the error satisfies 



(41) 



(42) 



(43) 



i<^<fc 

(T = 0,1 



Holder's inequality, applied to the first sum, gives 



1/2 



iv +1 . 



l<f<fc 



+ 1)' 



£=1 



and by (g^), (|3|), and (0), we finally get 

e{SN,AJ) < c{ky^n,^f\\L^ +2%' + n^H\\f\\lN + 1)) 
< crT^ log^^^ n loglogn max(||/||^iv, 1). 

This implies the statement for p = 2. 



□ 



Remark. Since quantum algorithms are not linear, the statement of The- 
orem 1^ does not give any information on / E Lp of norm greater than one. 
Our proof, however, does. It shows that the algorithm developed for fixed 
1 < p < oo and n, iV G N has the property that for all / G 



e{SN,AJ)<c{ 



n-imax(||/||^,,l) 

log ' n loglogn max 
n-2(i-Vrtmax(||/||^,,l) 



if 2 < p < oo 
if 1< p < 2. 
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4 Lower Bounds 



In this section we derive lower bounds on the quantities en{S, F) first in the 
general setting and then for F = Bp , S = Sn. Let D and K be nonempty 
sets, let L G N and let to each u = (uq, • • • , ul~i) G {0, 1}^ an fu G J^{D, K) 
be assigned such that the following is satisfied: 

Condition (I): For each t £ D there is an < £ < L — 1, such that 
fu{t) depends only on U£, in other words, for u, u' G {0, 1}^, ui = ti^ implies 

fu{t) = fu'{t). 

This type of function system will play a key role in our lower bound proofs. 
Condition (I) is easily seen to be equivalent to the following 

Condition (la): There are functions g^^gi G J-{D,K) and a decompo- 
sition D = U^^Q Di with De n A/ = (/ / /') such that for t G A 

f (i\ - j 9o{t) if Ui = 
I 9i{t) if ue = l. 

The first result is based on the polynomial method by Beals, Buhrman, 
Cleve, and Mosca (1998) and extends their Lemma 4.1 to our general setting. 

Lemma 9. Let L G N, let {fu)u&{o,i}^ — ^i^i^) ^ system of func- 
tions satisfying condition (I), and let A he a quantum algorithm on J^{D, K) 
without measurement, m = mq(A), n = nq(A). Then for all x,b £ Z[0, 2™), 
Af^{x, b) (defined in ^ and considered as a function ofu, is a complex 
multilinear polynomial in the variables uq, . . . ,ul-i of degree at most n. 

Proof. Let 

A = {Q, {U,)]^o), Q = (m, m',m", Z, r, /5). 
Fix b G Z[0, 2™) and define Wj and Pj{x, u) for j = 0, . . . , n by 

Wj = UjQuUj-iQf^ . . . UiQf^Uob = Pj(^^ ^) 1^) ■ 

xez[o,2™) 

Then 

Pnix,u) =Af^{x,b). (44) 
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Of course, po{x,u) are constants, so polynomials of degree in u. Now we 
proceed by induction over j. Assume that for some j, < j < n, the Pj{x, u) 
are polynomials of degree < j in u. Define qj{x,u) by 

QfuWj= Yl <lj{x,u)\x). 
xez[o,2'") 

Since 

Qfu'Wj = Qu Pjix,u) \x) = Y Pjix,u) Qf^ \x) , 

xGZ[0,2™) zeZ[0,2'") 

and since Qf^ is a bijection on the basis states, we get 

qj{x,u) =pj{Qj^x,u). 

Now fix X G Z[0,2"). Represent |x) as \i) \y) \z) with i G Z[0,2™'), y G 
Z[0,2™") and z E Z[0, 2™-"^'-™"). According to the query definition (|), 
we have Qj^ \x) = \i) \y) \z) if i ^ Z. Hence, in this case qj{x,u) = pj{x,u), 
so deggj(2;, ■) < j. U i £ Z, 

Q]^'\x) = \i)\yeP{fu(.rm)\z). 

Let, according to condition (I) above, i be such that < i < L — 1 and 
fu{T{i)) depends only on U£. We denote fu{T{i)) = so for U£ = and 
fu{T{i)) = si for Ui = 1. It follows that 

Q7i|x) = |0|ye/3(A(r«)))k) = | ^ = ? 

' ' ' ' '--\JU\ \ J/J/ I / y \t)\y Q p(si))\z) := Xi if Ue = l. 

Consequently, 

qj{x,u) = pj{Q^^x,u) = {1 - Ui)pj{xo,u) + UiPj{xi,u), 
which implies deg qj{x, • ) < j + 1. Now 

Wj+i = Uj+iQf^Wj = Uj+i Y 1j(y^ u)\y) , 

y6Z[0,2™) 

which gives 

Pj+i{x,u) = Y Uj+i{x,y)qj{y,u), 
j/ez[0,2™) 
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where {Uj+i{x,y))x^y^z[o,2-^) is the matrix of the transformation J/j+i in 
the canonical basis. Since the Uj+i{x,y) are scalars not depending on u, 
and since deg^j(x, • ) < j + 1, it follows that degpj+i(2;, • ) < j + 1. This 
completes the induction and shows that degpn{x, ■) < n. Now the lemma 
follows from (^) and the observation that, since the Uj take only the values 
and 1, we can replace any polynomial by a multilinear one without changing 
its values on {0, 1}^. □ 

Corollary 2. Let L G N and assume that (/n)ne{o,i}^ — ^{^^^) satisfies 
condition (I). Let A he a quantum algorithm from J^(D,K) to a normed 
space G. Then for each subset C G, 

p{u) = A{fu){G} 

is a real multilinear polynomial of degree at most 2nq[A). 

Proof. This follows readily from Lemma |^ and relations (j^ and (|9|) . □ 

The next lemma is based on the results of Nayak and Wu (1999). To 
state it, we introduce some further notation. Define the function £)(L, i') 
for L G N, < £ / f < L by 

Note that j{L - j) = {L/2f - (L/2 - jf, so this expression is minimized 
iff \L/2 — i\ is maximized. For u G {0, 1}^ set \u\ = "YldZ^ "^i- 

Lemma 10. There is a constant cq > such that the following holds: Let 
D,K be nonempty sets, let F C J^{D, K) be a set of functions, G a normed 
space, S : F ^ G a function, and L G N. Suppose ifu)ue{o,i}'- ^ ^{D,K) 
is a system of functions satisfying condition (I). Let finally < £ ^ £' < L 
and assume that 

fu^F whenever \u\ £{£,£'}. (46) 

Then 

eUS,F) >^mm{\\S{fu)-S{fu')\\ \ \u\ = £, \u'\ = i'} (47) 
for all n with 

n<cog{L,i,i'). (48) 
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Proof. Nayak and Wu (1999, Theorem 1.1) showed that there is a constant 
c > such that for all L G N and < i £' < L the following holds: If p is 
an L-variate real polynomial such that 

-1/4 < p(n) < 5/4 for all nG{0, 1}^, 

3/4 < j5(u) < 5/4 if u € {0,1}^, \u\=i, 

and 

-1/4 <p('u) < 1/4 if u G {0, 1}^, |n| 

then 

degp>cg{L,£,i'), (49) 



where g was defined in (45). Denote for j = £,£' 

G, = {SiU)\\u\=j} (50) 

and 

6 = d{Ge,Ge'), (51) 

where for X,Y <^ G, 

d(X,Y)= inf \\x-y\\. 

(For X £ G we write d{x,G) instead of d{{x},G).) Now let A be any 
quantum algorithm from F to G with nq{A) = n and 

e{S,A,F) <6/2. (52) 

As we mentioned after the definition, a quantum algorithm on F is always 
also a quantum algorithm on T{D,K). For each u G {0,1}^, let C,u be a 
random variable with distribution A[fy). Define 

p{u) = A{fu){g G G I d{g, Gi) < 6/2} = P{d(Cn, G,) < 5/2}. 

It follows that 

< piu) < 1 (53) 
and, by Corollary |^, p is a real polynomial satisfying 

degp < 2n. (54) 
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Because of (Mffl and ([5^), we have for |u| = ^, 



3/4 < P{\\Sifu)-Cu\\<S/2} 

< P{d{Cu,Ge)<5/2}=p{u). (55) 



On the other hand, for \u\ 



1/4 > P{||5(/„)-C«|| > V2} 

> P{d{Cu,Ge,)>S/2} 

> P{d{Cu,G,)<5/2}=p{u). (56) 



From (^ - 5f:) and (^), we infer 

2n > degp > cg{L, £,£'). 

Now choose any cq < c/2. Then n < coq{L,£,£') imphes en{S,F) > 5/2, 
which, because of (|50|) and (|5l|), is the same as (^7|). 

□ 

The following theorem is the main result of this section. The case p = 
oo is due to Nayak and Wu (1999), and the case 2 < p < oo is a direct 
consequence. For the sake of completeness we include this part in the proof 
below. (Another reason for this is that we use a slightly more general notion 
of query, so this way we formally check that their bound holds true also for 
our model.) 

Theorem 2. Let 1 < p < oo. Then there are constants cq,ci,C2 > such 
that for n,N eN, 

aro t^N^-^ f n-2(i-i/p) if l<p<2 and n < cqVN 
^ In if2<p<oo and n < ciN. 



Proof. Let cq be the constant from Lemma 10. Let 1 < p < 2 and 

n < cqVn. (57) 

Define 

L = [cq 2^2] , £ = 0, £' = 1. 
It follows from (|57| ) that 1 < L < A^. Moreover, 

n<coVL = coQiL,£,£') (58) 
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and 

L < c^'^r? + 1 < (cp ^ + \)r?. (59) 
Put M = \ L-^N\. Hence 1 < M < iV and 

M <L-^N < 2M. (60) 
Define if^j {j = 0, . . . , L — 1) by 



(iV/M)i/P if jM <i< {j + 1)M 
otherwise. 



V'j(i) = 

Note that ^pj G and 

For each u = {uq, . . . , ul-i) G {0, 1}^ define 

L-l 

fu = ^Uj'4}j. (61) 
i=o 

Since the functions V'j have disjoint supports, the system (/u)u£|o,i}^ satis- 
fies condition (I). Lemma 10 and relation (^) together with (|60| ) and (|59| ) 
give 

el^{SN,Bp) > ^ min {IS'Ar/u - 5Ar/„/| I |n| = 0, = 1} 

= ^ (MiV-i)'"'/^ > i(2L)-(i-i/p) > c2n-2(i-Vp) 

for some constant C2 > 0. This proves the statement in the first case. 

Now we consider the case 2 < p < oo. Since C Bp whenever p < oo, 
it suffices to prove the lower bound for p = oo. We set ci = 2-^(co ^ + 2)"^ 
and assume n < ciN. Let 

L = 2[co^n + l], £ = L/2-l, £' = £ + 1 = L/2. 

It follows that L > 4 and 



q{L, e, £') > min j{L - j) = ^L^/A - 1 > Cq ^n. (62) 

3=1,1' 
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Moreover, since 1 < n < ciiV, we get 

L = 2 \cQ^n + l] < 2(co + 2) < 2{cq^ + 2)n < N. (63) 
Now let M = [L-^N\ , then (|6|) holds again. Set 



1 if jM <i<{j + l)M 
otherwise 



for j = 0, . . . , L — 1, and let /„ be again defined by (^T|). Clearly, {fu)ue{o,i}'- 
satisfies condition (I) and /„ G for all u G {0, 1}^. Lemma |l^ together 
with relations (|6^), (|60|) and (|^) gives 



e^(S^,i3^) > ^mm{\SNfu-SNfu'\\\n\=£,\u'\=£ + l} 

= -MiV-i > — > C2n-i 
2 - 4L - 

for some C2 > 0. □ 

Remark. Comparing Theorem ^ with Theorem |l|, we see that matching 
upper and lower bounds were obtained except for the case of 1 < p < 2, 
n > cq^/N . This case is settled in Heinrich and Novak (2001b). 



5 Integration in Lp([0, 1]^) 

Here we present an application of the summation results to integration of 
functions. Further results will be contained in Heinrich (2001). Let 1 < p < 
oo, d G N, Z? = [0, 1]*^ and let Lp{D) denote the usual space of p-integrable 
with respect to the Lebesgue measure functions on D, equipped with the 
norm 

\\f\KiD) = ^l^\f{t)\'dt 
if p < cxD and 

ll/IUoc(D) = esssuptg£,|/(t)|. 
Let Id : Lp{D) ^ R be the integration operator, defined for / G Lp{D) by 

Idf= I f{t)dt. 

JD 
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In this chapter we will consider G = R and S = Id- We want to integrate 
functions from the unit ball B{Lp{D)) in the quantum model of computa- 
tion developed in section 2. Strictly speaking, Lp[D) consists of equivalence 
classes of functions being equal almost everywhere. Hence, function val- 
ues are not well-defined, in general. This changes, however, if we consider 
subsets of Lp{D) which consist of continuous functions, or more precisely, 
of equivalence classes which contain a (unique) continuous function. This 
is how we shall approach the integration problem - we study it for certain 
subsets 6 C B{Lp{D)). We shall assume that S is an equicontinuous set 
of functions on D. Since D is compact, equicontinuity is equivalent to uni- 
form equicontinuity, and the latter means that for each e > there is a 
6 > such that for s,t & D, \\s — i||oo < ^ implies \ f{s) — f{t)\ < e for 
all f E F. Note also that it follows readily from the Arzela-Ascoli theorem 
that £ C B{Lp{D)) is equicontinuous iff £ is relatively compact in the space 
C{D) of continuous functions on D, equipped with the sup-norm. (A simi- 
lar approach was chosen in Novak, 1988, to discuss restricted Monte Carlo 
methods.) 

Theorem 3. Let 1 <p < oo. Then there are constants ci,C2 > such that 
for all d, n G N 

cin~^ < sup e'^{Id,£) < C2n~^ 2 < p < oo 

£cB{Lj,{D)) 

cin"^ < sup e^{Id,£) < C2n~''^ log^/^ nloglogn 

£cBiL2{D)) 

cin-2(i-Vp) < sup el{Id,£) < C2n-2(i-Vp) l<p<2, 

£cB{Lp{D)) 

where the supremum is taken over all equicontinuous subsets £ of B{Lp{D)). 

Proof. First we prove the upper bounds. Let £ C B{Lp{D)) be equicontin- 
uous and let n G N. For A; G N let 

D= [j Di 

i=0 

be the partition of D into 2'^^ congruent cubes of disjoint interior. Let Sj be 
the point in Di with the smallest Euclidean norm. Let Pk be the operator 
of piecewise constant interpolation with respect to the partition (£'i)^=o~^ 
in the points (si)i=o~^ (to avoid ambiguity, if a point belongs to more than 
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one of the sets Di, we assign to it the value f{si) for the smallest such i). 
Due to the equicontinuity of £ there is a A; G N such that 

\\f-Pkf\\L^iD)<n-' (64) 

for all f ££. Fix this k and put N = 2'^''. It follows that 

suv\Idf -Id{Pkf)\<n-\ (65) 

Moreover, defining 

T:£^L^ by r(/)(i) = /(s,) (i = 0, . . . , iV - 1), 

we get 

N~l 

N 



N~l 

h{Pkf) = T7 E /(^*) = ° r(/). (66) 



i=0 

Note that for / G f C B{Lp{D)) 



Af-l \ Vp 



< 



p(D) + 11/ - Pkf\\Lj,{D) 



< \\f\\LAD)+n-' <2. 

Consequently, F maps £ into 2;S^. Lemma |6|, Corollary |^ and relations (|6 
and (|66|) imply 

(/,,£:) < n-i + e^(5^oF,£:) <n-i + e^(5iv,2^^) 
= n-i + 2e^(5;v,<), 

hence Theorem |l] yields the upper bounds. 

To verify the lower bounds, fix a o" with < a < 1 and let ip he a 
continuous function on R"^ with 

supp V Q [0, l]'^, < < 1 and /^V = o". 

Fix n E N and choose N = 2'^^ in such a way that 

cqVn > n and ciN > n, (67) 
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where cq and ci are the constants from Theorem Set 

^i{t) = ^{2\t - Si)) (i = 0,...,Af-l), 
with the Si as in the preceding part of the proof. Consequently, 

/rfVi = 2-*/dV' = f^2-* = aiV-i. (68) 

Define 

N-l 

T:B^^Lp{D) by T{f) = f{€)^i. 

For / G 



j=0 



N-l N-l 



N-l „ N-l 

= 2-'^^ ^ |/(i)|W |^(t)Rt<iV-i j;|/(i)|^<l. 



i=0 ^ j=0 

We define £ = T{Bp), which is a subset of B{Lp{D)). Since the functions 
V'i are continuous, and < N^^'p for all / G ^B^, the equicontinuity of £ 

easily follows. Furthermore, 

Af-l N-l N-l 

h o r(/) = hY. /(^)^* = E f^^^d^^ = E ^(^) = ^-^^^^ 

1=0 i=0 i=0 

Lemma ^ and Corollary |l| give 

ael{SN,B^) = eU^SN,B^) < el{I,,£), 
and the result follows from relation ( pTj ) and Theorem |2[ □ 

6 Comments 

Our results were formulated in the language of information-based complex- 
ity theory - the minimal error at given cost (number of function values, 
functionals etc., in our case queries). Lower bounds in terms of the number 
of queries mean the more that no algorithm can have better arithmetic (bit) 
cost. On the other hand, if we have upper bounds on the number of queries, 
this does not necessarily mean a corresponding estimate of the cost in the 
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bit model. However, for the problems considered in this paper we encounter 
a situation which is largely parallel to the experience in information-based 
complexity: As a rule, the developed algorithms, which are optimal in the 
query sense, show a similar behaviour (usually up to certain logarithmic 
terms) in their arithmetic (bit) cost. Let us have a closer look at our algo- 
rithms from this point of view. 

The bit cost of one query of the type (|) we define to be m' + m" (the 
number of bits to be processed). When we consider the bit cost, let us assume 
that both and n are powers of two, which is no loss of generality since the 
other cases can be reduced to that. We also assume n < N, see the remarks 
before Theorem [l|. The algorithm from Lemma makes one measurement 
and can be implemented on O(logA^) qubits using 0{n logA^) quantum 
gates. The algorithm of Lemma ^ requires 0{logN) qubits, 0{i'nlogN) 
gates and makes 0{i') measurements. Finally, the algorithm from Theorem 
m needs O(logA^) qubits, 0{nlogN) gates and 0(logn log log n) measure- 
ments for p < oo (one measurement if p = oo). 

To discuss the algorithm of Theorem ^, let us introduce the following 
quantity for an equicontinuous subset £ C B{Lp{D)) and e > 0: 

K{£,e) = minj/c G N | 

\f{s)- f{t)\ < e whenever / G .5, s, t e ||s - t||oo < 2"''}. 

Then for a given £ C B{Lp{D)) we have to compute the mean of = 
2dk jiujnbers, where it suffices to take k = K{£,l/n). If N < n, this 
can be done with 0{logn) qubits, 0{N logn) gates and one measurement 
(see the remarks before Theorem If n < A^, this can be implemented 
on 0{dK{£, 1/n)) qubits, with 0{dnK{£, 1/n)) gates and ©(logn log logn) 
measurements for p < oo and one measurement for p = oo. (The constants 
in the O-notation do not depend on £ and d.) 

Next let us compare the results obtained above to the classical deter- 
ministic and Monte Carlo setting. We denote the respective quantities by 
e'lf^ and e™'^. This discussion is carried out in greater detail in Heinrich 
and Novak (2001a), where also the related definitions and references can 
be found. The following table contains the order of the respective quanti- 
ties, that is, the behaviour up to constants. We also omitted the additional 
logarithmic factor in the case p = 2. Furthermore, we assume for the case 
Bp that n < ciN, where in the classical settings, ci is any constant with 
< ci < 1, while in the quantum setting for 2 < p < oo, ci is the constant 
from Theorem |2[ Moreover in the quantum setting for 1 < p < 2, we assume 
n < co^/N, with cq from Theorem ^, as well. Finally, when we write Bl^, we 
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mean (in all three setttings) the supremum over all equicontinuous subsets 
S C B(Lp([0, l]*^)) as in the previous section. 





pdet 








2<p<oo 


1 




n" 


-1 




l<p<2 


1 




n~ 


-2+2/p 




2 < p < cx) 


1 




n~ 


-1 




l<p<2 


1 




n~ 


-2+2/p 



The result on Bl^ in the randomized setting can be found in Heinrich (1993). 
The respective statement for the deterministic setting is easily derived using 
standard methods of information-based complexity theory. A little further 
below we indicate the proof of a somewhat stronger result. 

It might be illustrative to formulate the results in terms of complexity. 
Here we impose the corresponding restrictions. We always assume e < eo 
for some constant > 0. In the quantum setting, the case 1 < p < 2 holds 
only for > c(l/e)^/*-^~^\ for some constant c > 0. Again, the case p = 2 
holds up to logarithmic terms. 





compi^ 


comp^'^ 


compi 


B^,2<p<oo 
B^,\<p<2 

Hl^, 2<p< DO 

Bl,,1<p<2 


A^ 
A^ 

GO 
GO 


min((l/£)2,A^) 
min((l/£)*'/(P-^,A^) 
[l/ef 
(l/e)?>/(p-i) 


min((l/£),A^) 
min((l/£)?'/(2(?'-i)),A^) 

(l/e)p/(2{p-i)) 



In the case Bl-^ we have oo in all three settings. For Bi we have A^ in both 
classical settings, while in the quantum setting our results give the lower 
bound \/N and the (trivial) upper bound N. The question of the correct 
order of complexity in this case is answered in Heinrich and Novak (2001b). 

We see that for the problems considered here quantum algorithms reach 
a quadratic speedup over classical randomized ones and - at least as far 
as the pure number of queries is concerned (disregarding the bit cost and 
number of qubits) - an arbitrarily large speedup over classical deterministic 
algorithms. Let us discuss this last point in some more detail and also 
address the bit issue again. Namely, we show that there are equicontinuous 
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sets £ in ;B(Loo([0, 1])) with arbitrarily slowly decreasing e'^^{I\,£). More 
precisely, for any sequence (£n)n£N with 



< e„ < 1, En+i <en< 2e2n, and lim e„ = (69) 

n^oo 

there is an equicontinuous set £ C ;B(Loo([0, 1])) such that for all n G N 

ei'\h,£)>en/?>2. (70) 

Indeed, we define £ as the set of functions / on [0, 1] such that for all G N 
and s,t E [0,1], |s -t| < implies |/(s) - f{t)\ < £2^- Let 



t if < t < 1/2 

(1-t) if l/2<t<l 
otherwise, 



e^f-i(/i,f) > £2^/8 > £2^-1/16, 



and put for G N and < « < 2^^ — 1 

i^kAt) = e2^^{2\t - 2~h)). 
It is easily checked that for any a, G {—1, 1} (^ = 0) • • • > 2^^ — 1), 

2*-l 

i=0 

and Iiipk,i = 2~('^+^)e2fe- ^ standard argument from the deterministic set- 
ting of information-based complexity theory (see e.g. Novak, 1988, Prop. 
1.3.5 b) yields 

-2^ 

which implies ([70|). Recall, on the other hand, that by Theorem ^, 
ei{Ii,£) < cn~^. 

Now let us turn to the bit cost. We show that an exponential speedup 
is possible. Fix any 7 with < 7 < 1. We choose £1 = 1 and £„ = 
(logn)"''' (n > 1). This sequence satisfies (|69[). Let £ C B{L^{\}},1])) be 
the corresponding set constructed above, so that 

ef*(/i,£:)>(logn)-732 (n > 1), 

which means that for any e with < £ < 1/32 we need at least 2^"^/*^^^'^'')^''^ , 
that is, exponentially many operations to reach error e deterministically. By 
the construction of the set £ we have 

K{£,l/n) < [n^/^l, 
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which imphes, by the discussion at the beginning of this section, that in 
the quantum setting, an error of e can be reached with 0{\/e) queries, 
one measurement, 0{{\/ e)^/'^) qubits and 0{{\/ eY^'^^^) gates, that is, with 
polynomial total cost. 

Finally we discuss a topic concerning the relations to information-based 
complexity. A look at our notion of a query might lead to the impression that 
it covers only what is called standard information, that is, function values of 
/, while in information-based complexity also more general types of informa- 
tion are considered (e.g. arbitrary linear functionals or scalar products with 
certain basis functions). This could be relevant not only in finite element 
methods, but also in the case that function values are not well-defined. Let 
us show how our approach covers also this situation. 

So let F and K be nonempty sets, 5 : F — > G be a mapping from F 
to a normed space G and let A be a nonempty set of mappings from F to 
K. We seek to approximate S again, but now the algorithm is supposed to 
use information about f G F oi the form A(/) for A G A. Let us define a 
A-based quantum algorithm from F to G to be simply a quantum algorithm 
A from J^(A, K) to G. Introduce the mapping 

^ -.F ^ :r(A, K) 

defined for f £ F hy 

^(/)(A) = A(/) (AG A). 

The error of ^ at / G F is defined as follows. Let be a random variable 
with distribution A{^{f)). Put 

e{S,A,f,0)=mi{e \ P{||5(/) -C|| > e} < 6} . 

Various further quantities like e{S,A,F), en{S,F,A) etc. can be defined on 
this basis as in section 2. The results of section 2 as well as the general results 
of section 4 remain valid for this situation if formulated appropriately, that 
is, if applied to A as an algorithm from J-'{A, K) to G. The resulting form of 
the unitary mappings associated with the query is worth mentioning: Let Q 
be one of the queries being part of A. Since A is an algorithm on ^(A, K)^ 
its queries have the form ([l|), where everything is as specified there except 
that 

r : Z ^ A. 

Let us denote Aj = T{i) for i £ Z. Then an element f G F gives rise to the 
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following unitary operator implementing the query 

n u\w\u,\-\ N)l^©/5(A.(/)))|y> if iG^ 

\i) \x) \y) - I I otherwise. 
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